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Abstract 

In  [1]  (resp.  [2])  we  studied  the  pollution-error  in  the  /i-version  of  the  finite 
element  method  and  its  effect  on  the  quality  of  the  local  error  indicators  (resp.  the 
quality  of  the  derivatives  recovered  by  local  postprocessing)  in  the  interior  of  the 
mesh.  Here  we  show  that  it  is  possible  to  construct  a-posteriori  estimates  of  the 
pollution-error  in  a  patch  of  elements  by  employing  the  local  error  indicators  over 
the  entire  mesh.  We  also  give  an  adaptive  algorithm  for  the  local  control  of  the 
pollution-error  in  a  patch  of  elements  of  interest. 
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1  Introduction 


Let  ft  C  denote  a  polygonal  domain  with  boundary  5ft  =  F  =  U  Fi, 

>=1 

where  F,,  i  =  1, . . . ,  Af  are  the  open  straight  edges  connecting  the  endpoints  Ai 
and  i4,+i  (Ai  =  (see  Fig.  1).  We  will  denote  the  internal  angles  at  the 

vertices  by  y?!, . . . ,  (0  <  <  27r;  if  (p,  =  27r  we  have  a  slit-domain  and  in 

this  case  we  have  to  understand  the  boundary  in  the  obviously  modified  sense). 
Further,  let  F  =  Fd  u  F;v,  f^)  n  Fyv  =  0  where  F/j  is  the  Dirichlet-  and  Fyv  is  the 
A^ewmann-boundary.  We  will  consider  the  mixed  boundary-value  problem  for  the 
Laplacian: 


—An  ;= 

—V • Vu  =  0 

in  ft 

(la) 

u  =  0 

on  Fd 

(16) 

du 

dn 

:=  Vu  -  n  =  g 

on  Fat 

(Ic) 

Here  u  is  the  exact  solution;  n  denotes  the  exterior  unit-normal  on  Fat  ;  5  €  L‘^{Tn) 
and  is  smooth  in  each  Fj ,  j  =  1, . . . ,  M. 

The  variational  formulation  of  the  model  problem  (1)  is; 

Find  u  €  /fpc  •—  |  u  =  0  on  Fcj  such  that 

Ba(u,v):=  f  Vu-Vv=  f  gv  V  u  e  (2) 

If  Fp  =0  it  is  assumed  that  g  satisfies  the  consistency  condition  J  p  =  0;  in  this 

case  the  solution  u  is  determined  uniquely  up  to  an  arbitrary  constant. 

For  the  finite  element  method  we  partition  the  domain  ft  into  triangular  or 
square  elements  with  straight  edges  r  defined  by  the  mesh  and  let 

s;(n)  e /ff J  €  5J(t)  Vren}  (3) 

Here  S^(r)  denotes  the  finite-element  space  over  t  and  p  is  the  degree  of  the 
elements.  For  the  meshes  of  triangles  5^(t)  =  ‘Pp('r),  while  for  the  meshes  of 
squares  S^(t)  denotes  the  tensor-product  polynomial  space  of  degree  p 

Shi'r)  I  ^2  }  (4) 

0<«iJ  <P 

The  finite  element  approximation  of  the  solution  of  (1)  satisfies: 
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(5) 


Find  Uf^  €  :=  fl 

Bq{u^,v,,)  =  f  gvk  y  v^e  5j,r^ 

There  exists  a  large  literature  devoted  to  the  a-posteriori  estimation  of  the 
error  of  the  finite  element  solution  ef^  :=  u  —  u^,  and  to  the  recovery  of  the  data  of 
interest  (e.g.  derivatives,  stresses,  etc.).  For  a  complete  discussion  and  references 
to  the  vast  literature  see  [3]- [8]  for  the  problem  of  a-posteriori  error  estimation  and 
[9]-[13]  for  the  problem  of  the  recovery.  All  of  these  procedures  have  local  character 
i.e.  the  computed  data  in  the  element  or  patch  of  elements  depend  on  the  finite 
element  solution  and  the  input-data  only  in  the  neighborhood  of  the  element  or 
patch  under  consideration.  In  [6]- [8]  and  the  theory  presented  there  we  underlined 
the  assumption  that  the  mesh  outside  the  region  of  interest  has  to  be  appropriate 
so  that  there  is  no  significant  pollution,  which  cannot  be  detected  by  any  local 
analysis.  We  have  shown  in  detail  the  effect  of  the  pollution  on  the  a-posteriori 
error  estimation  and  recovery  in  [1]  and  [2],  respectively. 

In  [1]  and  [2]  we  have  seen  that  the  pollution-effect  is  negligible  for  meshes 
which  are  constructed  adaptively  with  respect  to  the  energy  error-measure  in  the 
entire  domain  (we  will  call  such  a  mesh  a  fully- adaptive  mesh).  Nevertheless,  very 
often  in  engineering  computations,  the  mesh  employed  is  very  different  from  a 
fully-adaptive  mesh  because  of  limited  computer  resources.  It  is  often  assumed 
that  the  computed  data  are  reliable  (if  this  is  indicated  by  the  local  error  estimates 
or  the  difference  between  the  recovered  data  and  the  data  computed- directly  from 
the  finite  element  solution)  without  paying  attention  to  the  refinement  of  the  mesh 
outside  the  region  of  interest,  as  it  is  done,  for  example,  in  the  global- local  approach 
[14]-[16].  This  assumption  can  be  correct  or  false,  as  has  been  shown  in  the  results 
given  in  [1]  and  [2]  (see  also  [17]).  From  these  results  it  is  clear  that  there  is 
urgent  need  for  the  a-posteriori  control  of  the  pollution.  This  control  cannot  be 
local,  has  to  be  cheap  and  accurate  and  will  be  addressed  in  the  present  paper. 
We  will  assume  that  we  have  at  our  disposition  error  indicators  which  are  very 
robust  modulo  the  pollution  (i.e.  the  element  effectivity  indices  are  close  to  one 
if  there  is  no  pollution).  We  have  shown  (see  [6]- [8])  several  examples  of  robust 
indicators  in  the  interior  of  the  mesh.  The  analysis  of  the  various  error  indicators 
for  the  elements  at  the  boundary  and  especially  for  the  unsmooth  solutions  (in  the 
neighborhood  of  the  singular  points)  will  be  done  in  forthcoming  papers.  In  the 
present  paper  we  mention  only  a  few  illustrative  examples  related  to  the  standard 
estimators. 


2  Local  error  indicators 

Below  we  describe  two  types  of  local  error  indicators  which  will  be  used  in  the 
numerical  examples  for  the  construction  of  a-posteriori  estimates  of  the  pollution. 
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2.1  Error  indicators  based  on  patch-recovery 

Let  ujx  denote  the  patch  of  elements  connected  to  vertex  X.  For  each  patch 
Ux  we  recover  the  smoothened-gradient  ,  by  solving  the  following  least-squares 
problem; 

where  denotes  a  set  of  sampling  points  in  u;;^  (given  in  [18-20])  and 

nip  r  2  ■ 

S  E{9i(ynx))  (7) 

m=l  L»=l 

A  continuous  piecewise  polynomial  recovered  gradient  is  obtained  over  the 
mesh  by  combining  the  smoothened-gradient  q^  from  the  patches  which  correspond 
to  the  vertices  of  the  element  (see  [18-20]  for  the  details).  In  the  construction  of 
q^^  we  employ  only  patches  uix  for  vertices  X  in  the  interior  of  the  domain.  The 
values  of  q^^  at  the  boundary  are  obtained  using  extrapolation  from  the  interior 
patches.  The  element  error  indicators  are  given  by: 

r]r^  '■=  \\q^^ (8) 

In  the  numerical  examples  we  will  call  this  estimator  the  ZZ  (Zienkiewicz-Zhu)- 
estimator. 


2.2  Error  indicators  based  on  element-residual  problems 

Let  e  be  an  edge  and  denote  the  jump  of  the  normal  derivative  ofuf^  on  e, 
defined  by: 


Here  n,  and  denote  the  unit-normal  and  the  elements  associated  with  the 
edge  e,  as  shown  in  Fig.  2.  We  will  denote  by  r^  :=  Au^|  the  interior  residual  for 
element  t  and  let  — »  R  be  the  residual  functional  given  by 


D 


(10) 
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The  error  e;^  satisfies  the  residual- equation: 

Find  6;,  6  such  that 

Bn{e,,v)  =  •£  ^Av)  (11) 

'6T. 

Let  us  assume  that  the  element-residuals  have  been  modified  in  the  following 
way  (for  the  construction  see  the  Appendix): 

+  Y.  f  <  (12) 

such  that 

=0  V  V  e  sr(T)  (13) 

Here  6‘  is  the  correction  for  the  edge  e  and  the  element  r  and  denotes  the  set 
of  interior  edges.  For  any  interior  edge  e  the  correction  0^  is  constructed  such  that 
9i  =  —9i  .  We  then  have 

'out 

Y,  :fm  =  E  (14) 

Let  us  introduce  the  local  problems: 

For  each  element  t  find  €  H^{t)  such  that 

=  Vt,£/f>(T)  (15) 

Note  that  exists  and  is  unique  up  to  a  constant  since  =  0.  (If  5Tr)r£)  ^ 

0  we  let  =  0  for  edges  e  €  ^rflFp  and,  in  this  case,  the  solution  is  unique.) 

Let  us  denote  by  |l|v|||5  :=  •\JBs{v,v)  the  energy  norm  over  any  5  C  fi.  We  then 
have 


lllefcllln  =  sup 


:£h _ 

llkllln 


sup 


E  :Fr{v) 

£h _ 

llbllln 


sup 

v€i/i 

D 


^ _ 

IIHIIn 


(16a) 


and  hence 


'fcllln=  sup 
v^m 

•  n 


E  B,{e,,v) 

_ 

IIMII 


T€r, 


(166) 


6 


Thus  an  upper-estimator  for  the  global  energy-norm  of  the  error  is  obtained  pro¬ 
vided  that  the  local  problems  are  solved  exactly.  We  will  denote  this  estimator  by 
ERp. 

In  the  implementations  we  employed  approximate  solutions  of  the  local  prob¬ 
lems  which  are  computed  as  follows: 

Find  G  such  that 

V  u  6  (17) 

with  the  corresponding  element  error  indicators 

_  |||4r+‘>|||,  (18) 

In  the  numerical  examples  we  will  denote  the  estimator  defined  by  (17)  and  (18)  as 
ERpPp  +  k  (element-residual  with  p-order  equilibration  and  (p -f-  fc)-degree  polyno¬ 
mial  space).  For  further  details  about  element- residual  estimators  see  [8],  [21-29] 
and  the  references  therein. 

Remark  2.1:  By  letting  k  — >  oo  in  (17)  we  recover  the  exact  solution  of 

the  local  problem  (15).  In  most  cases,  by  letting  A:  =  3  we  obtain  element  error 

indicators  which  are  practically  the  same  as  the  indicators  corresponding  to  the 
exact  solutions  of  the  local  problems. 

We  also  considered  the  following  variation  of  the  local  residual  problem: 

Find  G  Mp’^{T)  such  that 

Vt;GMP+'(T)  (19) 

Here  M^'*’^(t)  is  the  bubble-space  of  degree  (p-|- 1)  (see  [28]-[29]).  In  the  computa¬ 
tions  we  employed  the  hierarchical  element  families  given  in  [30]  and  let 
be  the  span  of  the  edge  shape-functions  of  degree  (p  -|-  1)  only  (i.e.  three  (resp. 
four)  edge  shape-functions  for  a  triangular  (resp.  square)  element)  together  with 
<ill  interior  shape-functions  in  *S'^''’^(r).  We  will  call  the  estimator  based  on  (19) 
the  ERpB-estimator  (element-residual  with  p-order  equilibration  and  with  (p-f  1)- 
degree  bubble  space).  (For  a  detailed  study  of  the  element-residual  estimators  in 
the  interior  of  the  mesh  see  [6-8].) 

We  also  considered  the  following  patch-residual  estimator:  Let  r  be  the  element 
of  interest,  be  a  patch  of  elements  which  includes  r  and  a  few  mesh- layers  around 
it  and  consider  the  local  problem: 

Find  such  that 

V  u  e  sioM  (20) 
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Here  we  let 


5?.oK)  :=  I  t’l,  €  ^Hr),  V  =  0  on  5u;,  n  }  (21) 

The  error-indicator  in  the  element  r  is  computed  by 

=  lll«Lr‘’ijll.  (22) 

We  will  use  the  notation  PRpPp  +  A:  to  refer  to  the  patch -residual  problem  with 
equilibration  of  order  p  and  (p-t-A:)-degree  polynomial  space.  Below  we  will  consider 
patches  which  consist  of  a  few  mesh-layers  of  elements  around  the  singular  point. 


2.3  The  quality  of  the  error-indicators  in  the  neighbor¬ 
hood  of  the  singular  points 

We  now  give  sample  results  on  the  quality  of  the  error  indicators  in  the 
neighborhood  of  singular  points  for  uniform  meshes  of  triangles  and  squares;  we 
will  present  a  more  systematic  study  in  a  forthcoming  paper. 

Let  us  consider  the  problem  (1)  in  the  T-shaped  domain  (—1,1)^—  [0, 1]  x  [0,  —1] 
(resp.  the  slit  domain  (-1,1)^-  [0,1])  and  assume  that  the  boundary-conditions 
are  such  that  (for  a  complete  study  of  the  properties  of  solutions  in  polygonal 
domains,  see  [31]-[34])  u(r,t?)  =  r“^(at?)  where  (r,  i9)  are  polar  coordinates  with 
r  =  0  at  the  singular  point.  We  employed  homogeneous  boundary-conditions  on 
the  edges  of  the  domain  emanating  from  the  singular-point  at  (0,0)  such  that 
V’(at?)  =  sin(at?)  (resp.  =  cos(aj?))  for  homogeneous  Dirichlet  or  mixed 

(resp.  homogeneous  Neumann)  boundary  conditions,  and  Neumann  boundary- 

conditions  computed  from  the  exact  solution  in  the  rest  of  the  boundary.  We 

12.  .  / 
employed  the  values  q  =  —  and  —  in  the  case  of  the  T-shaped  domain  I  resp.  a  = 

3  3  \ 

-  and  a  =  -  for  the  slit-domain^ .  We  computed  the  finite  element  solution  using 

the  uniform-meshes  of  triangles  with  h  —  —  shown  in  Fig.  4a,  for  the  T-shaped 

In 

domain,  and  in  Fig.  5a  for  the  slit-domain. 

In  Tables  la,  lb  (resp.  Tables  Ic,  Id)  we  give  the  effectivity  indices  for  the  vari¬ 
ous  estimators  in  the  neighborhood  of  the  singular  point  for  linear  (resp.  quadratic) 
triangles.  The  effectivity  indices  are  given  for  the  elements  near  the  singularity 
shown  in  Figs.  4b,  5b.  In  Tables  la-ld  we  included  the  estimators  ERpPp  -t-  1, 
ERpB,  ZZ  aind  PRpPp  -I-  3  (for  the  patch- residual  estimator  we  employed  two 
mesh-layers  around  the  singular  point  to  define  the  patch).  In  Table  2  we  give  the 
effectivity  indices  for  the  estimators  ERpPp  -f  A:  for  A:  =  1,  2,  3,  4,  PRpPp  -t-  3  and 

ZZ  for  the  T-shaped  domain  with  a  =  -  where  the  finite  element  solution  was 
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computed  using  the  uniform  mesh  of  bilinear-squares  shown  in  Fig.  6a  (j)  =  1^  h  = 

— ^ .  The  effectivity  indices  are  reported  for  the  elements  in  the  neighborhood  of 

the  singular  point  shown  in  Fig.  6b.  Note  that  for  the  meshes  of  squares  we  em¬ 
ployed  three  mesh-layers  around  the  singular-point  eis  the  patch  for  the  estimator 
PRpPp  -f  3. 

We  see  that  the  error  in  the  layer  of  elements  adjacent  to  the  singularity  can  be 
relatively  large  (e.g.  more  than  50%  of  the  error  in  the  entire  mesh).  The  elements 
where  the  energy  of  the  error  is  more  than  10%  of  the  global  energy  of  the  error 
will  be  called  key  elements.  The  large  spread  of  the  effectivity  indices  in  the  key 
elements  is  in  contrast  to  the  robustness  of  the  estimators  for  smooth  solutions  and 
interior  elements  that  we  reported  in  [6]-(8]. 

As  we  said  above  the  quality  of  the  pollution  estimates  depends  on  the  quality 
of  the  estimator.  Hence  we  should  also  consider  another  estimator  which  is  very 
accurate  in  the  key  elements.  Such  an  estimator  can  be  obtained  by  assuming  the 
aisymptotic  expansion 

3 

u(r,?9)  =  sin(Q,i?)  (23) 

i=l 

(or  «(r,  t?)  =  ^  C,  r“*  cos(Q,  t?)  depending  on  the  boundary-conditions]  and  by 
extracting  the  coefficients  C,  from  the  finite  element  solution  by  employing  the 
contour-integral  method  given  in  [30].  For  example  for  the  case  of  the  strongest 

singularity  a  =  -,  the  mesh  shown  in  Fig.  5a  and  the  elements  shown  in  Fig.  5b 
the  effectivity  indices,  obtained  by  the  extraction  method,  are. 


..  =  0.91, 

=  0.91, 

=  0.93 

=0.95, 

=  0.97, 

=  0.99 

(p  =  1 )  and 

=  0.91, 

=  0.92, 

=  0.94 

=  0.97, 

=  0.99, 

=  1.00 

{24a) 


{24b) 


for  quadratic  triaingles  (p  =  2).  Comparing  the  above  effectivity  indices  with  the 
effectivity  indices  reported  in  Tables  lb.  Id  we  see  significant  improvement  in  the 
quality  of  the  estimator  in  the  key-elements. 

Remark  2.1.  In  the  examples  below  we  know  the  exact  solution  and  we  could  get 
the  exact  error  indicators.  We  employ  the  extraction  method  as  a  general  approach 
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for  the  error  computation.  We  will  discuss  the  construction  of  error-indicators  in 
the  neighborhood  of  singular  points  in  a  forthcoming  paper. 

3  A  posteriori  estimate  of  the  pollution  error 

In  Section  2  we  introduced  various  error  estimates.  We  will  now  consider  the 
estimator  ERp. 

Let  uJf^  be  a  patch  of  elements  (possibly  consisting  of  one  element)  and  let  W/j 
denote  a  patch  which  consists  of  u)f^  and  a  few  mesh-layers  around  it  (e.g.  Fig.  3). 
From  (11)- (14)  we  get 


=  E  (25) 

Let  for  j  =  1,  2,  be  defined  by 

=  E  V  ^  e  (26) 

’■£«h 

Ba(K-,v)  =  E  V  €  Wf  (27) 

Obviously 

Vf*  -f  =  e,  (28) 

In  [1]  we  showed  that,  in  general,  we  have 

IIIK^NIU  as  /^-.O  (29) 

The  symbol  «  means  that  there  exist  constants  C^'',  which  depend  on  the 
geometry  of  the  mesh  in  but  not  on  the  mesh-size  fi,  such  that 

(30) 

In  [6]- [8]  we  computed  the  values  of  for  several  estimators  and  a  variety 

of  meshes  and  identified  several  robust  estimators  (e.g.  estimators  with  C'^'  and 
close  to  on**  for  the  types  of  meshes  which  occur  in  practical  computations)  for 

mesh-cells  in  the  interior  of  the  mesh.  Here,  as  usual,  we  let  •=  (  * 

with  T]r  being  the  element  error-indicator  corresponding  to  the  particular  estimator 
(e.g.  ERp).  From  (28)-(30)  (and  the  results  in  [1])  we  see  that  |1|V"2'*^'‘||L  is  the 
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measure  of  the  pollution  and  that  the  estimate  of  the  pollution  is  equivalent  with 
the  estimate  of  |||V2^'‘1||^,  .  We  also  note  that  is  harmonic  in  W/,. 

h 

Let  *  €  n,  be  the  function  satisfying 


= 


88 


dX: 


(*) 


0 

0 


dn' 


in  n 

on 

on  r^v 


(31a) 

(316) 

(31c) 


88 


Here  i  =  1  or  2  and  -r — (sc)  denotes  the  x. -derivative  of  Dirac’s  delta  centered  at 
ax, 

X.  Eq.  (31a)  is  understood  in  the  sense  of  the  theory  of  distributions.  Obviously 
^  -  By  the  standard  arguments  we  see  that 


dv^ 

8xi 


(*)  =  E 


(32) 


’■2*h 


Let  us  now  replace  the  function  G,-*^  by  the  difference 


(33) 


where  n,  denotes  the  unit-vector  in  the  i-th  coordinate  direction  and  G^^^  is  the 
(classical)  Green’s  function.  Then  for  any  x  €  fi  —  W/,  we  have 


|g!'I(*)  -  Gi«(x)|  <  Ch  (34) 

We  note  that  we  can  replace  G;  by  another  combination  of  Green’s  functions  so 
that  the  error  will  be  of  order  h'  (s  >  0,  arbitrary)  instead  of  order  h  as  in  (34). 
Using  (13)  and  (32)  we  get 


^(*)  =  E  -  <-r)  (35) 

t€T„ 

’■**h 

where  is  the  best-approximation  of  G|*^  on  r  from  5^(t)  in  the  energy-norm. 
By  employing  the  definition  of  (see  (15))  we  see  that 

8V^^ 

=  T,  -  w,)  (36) 

r€T^ 
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Hence 


^(i)|<  E  lll^,IIUi|Gr-«>,lll,  (37) 

’■2"h 

Replace  now  by  and  by  is  the  best-approximation  of  G,-^^  on 

to  obtain 

|^(*)|<  E  7,(«»)ll|Gj'’-<i,ilUl+a)  (38) 

’•Z*h 

where  T)^{Uf^)  is  the  ERp-indicator  computed  from  Uf^.  Because  Wj.  is  the  best- 
approximation  on  T  of  the  function  G-*^  we  get 

lllGl*'  -  -irlll,  <  -!,(GS)  (39) 

and  using  (38)  we  obtain 

|^(*)|  <  (  E  IrM  Vr(&S))(l+Ch)  (40) 

’•*4 

Here  we  are  assuming  that  *  is  a  nodal  point  of  the  mesh  and  that  TfA^l^h)  is 

the  ERp-indicator  computed  from  the  finite  element  solution  G|^^  ^  which 

satisfies 


Ba{Gi%v)  =  ^ (u(*  -I-  n,- h)  -  v{x))  Vue  (41) 

We  assume  that  x  +  n^h  is  also  a  nodal  point  of  the  mesh  (this  can  be  always 
achieved  with  obvious  modifications  of  the  approximation  of  G,-*^). 

We  note  that  G\^iI  approximates  well  the  function  G-*^  on  fl—  u;^,  although  G,-*^ 
does  not  belong  to  .  Note  that  has  singular  behavior  at  the  corner-points. 
For  further  details  see  [35]- [36].  Note  also  that  Gj|^^  (and  hence  the  corresponding 

error  indicators  is  cheap  to  compute  because  it  requires  only  the  solution 

of  the  system  of  discrete  equations  for  an  additional  right-hand  side. 

We  employ  the  indicator  T]^{G\^h)  ^  estimate  of  the  error  in  the  best- 

approximation  on  the  left  hand  side  of  (39)  and  hence  the  possible  pollution  does 
not  play  any  role.  Also,  from  the  construction  of  ^2'^'’,  we  see  that  it  relates  only 
to  the  residuals  (which  are  local)  and  hence  the  absence  of  the  pollution-term 
from  does  not  influence  the  accuracy  of  the  estimate  (40).  Hence  the  only 
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requirement  for  the  reliability  of  (40)  is  the  accuracy  of  the  indicators  modulo 
the  pollution.  Further,  since  is  smooth  in  we  can  approximate  |||V2^''11L 

h 

assuming  that  constant  in  by  letting 

(42) 

with  ®  €  u;^. 

Remark  3.1.  Although  we  mentioned  the  pollution  error  in  connection  with  the 
error  estimate  in  W/,,  the  same  idea  can  also  be  used  to  estimate  the  effect  of  the 
pollution  on  the  error  in  recovered  quantities. 

In  summary  the  a-posteriori  estimate  of  ^  (x)  is  constructed  tis  follows: 

'  OXi 

We  compute  the  indicators  T|^{uf^),  define 

^lf'(*)=  Y:  |B,  («.(«»),  I  (43) 

’•2"h 

where  and  are  the  error-indicator  functions  corresponding  to 

and  respectively, 

Mf\x)  =  Y  IrM  4.(G'!S)  (44) 


A4p'(*)=  E  (’!.(“*))=  E(’f.(G!S))^ 


We  then  have 


<  M\^\x)  (46) 

We  will  define  the  pollution-indicators  fi^^{x),  m  =  1,  2,  by 

,-!,V(i)  :=  I B,  (e.(«j),  4.(Gi5))  I .  ■=  ’!r(“-)  1. (<?&’)  (47) 

We  then  have  Ai\^\x)  =  ^  for  i  =  1,  2,  m  =  1,  2.  We  will  also  define 

pollution-estimates  for  |Vl^‘^'‘|(*)  in  the  form 
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m  =  1, 2, 3 


(48) 


Al'"'!*)  :=  V'(A4S"'(*))2  +  (M"'(*))^ 

Remark  3.2.  Although  the  derivations  are  valid  only  for  the  case  where  the 
error-indicators  are  based  on  the  ERp  estimate  we  can  expect  that  the  above  es¬ 
timates  give  reliable  answers  if  they  are  computed  using  indicators  other  than 
ERp-indicators  (provided  that  these  indicators  are  reliable  modulo  the  pollution). 

4  The  quality  of  the  a-posteriori  estimates  of 
the  pollution-error 

We  will  give  sample  results  on  the  a-posteriori  estimation  of  the  pollution  error. 
We  employed  the  estimates  given  in  Section  3  in  the  following 

form: 


M 


(0 

k,  aing 


-1-  M 


(«) 

k,  smooth 


(49) 


By  (resp.  \mooth)  denote  the  contributions  to  from  the  ele¬ 

ments  adjacent  to  the  singularity  (resp.  the  elements  in  the  rest  of  the  mesh).  For 
the  meshes  of  triangles  (resp.  meshes  of  squares)  we  computed  from  the 

first  layer  (resp.  the  first  two  layers)  of  elements  around  the  singular  point.  In  the 
results  the  computed  using  one  of  the  indicators  given  in  Sec¬ 

tion  2  (e.g.  ERpPp-l-3,  ERpB,  ZZ)  while  the  vVffc'^„„^-part  was  computed  either  by 
employing  an  extraction  method  or  by  employing  the  saime  indicators  used  in  the 
interior.  We  also  computed  estimates  for  |VV^'‘(*)|  by  combining  the  estimates 
for  the  ij-  and  Xj-components. 


dVn'' 

In  the  numerical  results  we  also  computed  the  “exact”  value  of  -~—{x)  by 

ax, 

solving  the  residual  equation  (26)  for  (a;)  using  elements  of  degree  (p  -f-  2) 


and  by  letting  (see  also  [1]) 


dx; 


avr- 

-gJ-(*)  _  — (»)  -  -gJ-(»,) 


(50) 


We  measured  the  quality  of  the  estimators  by  computing  the  effectivity  indices 

Al(">  ...  A(l") 


(m) 

k)  = 


i  =  1,2, 


:= 


dx, 


-(*) 


|V^"'‘(*)| 


(51) 


for  m  =  1,  2,  3. 
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In  Table  3  we  give  an  example  of  how  the  quality  of  the  estimates  of  the 
pollution  is  affected  by  the  quality  of  the  error-indicators  near  the  singular  point. 

We  considered  the  T-shaped  domain  with  ~  computed  the  pollution 

estimates  at  *  =  (0.5, 0.5)  for  the  uniform  mesh  of  triangles  shown  in  Fig.  5a  (h  = 

— )  for  linear  and  quadratic  triangles.  We  computed  the  pollution-estimates  by 
16' 

employing  the  error-indicators  from  the  ERpB  estimator.  The  results  in  Table 
3  show  overestimation  by  a  factor  of  1. 9-2.5  by  the  and  pollution- 

estimates.  This  is  due  to  the  fact  that  the  error-indicators  for  the  ERpB  estimator 
overestimate  the  error  in  the  key  elements  (see  the  results  in  Tables  la,  lb). 

In  order  to  be  able  to  evaluate  the  performance  of  the  pollution-estimates  in 
the  ideal  case  (without  the  influence  of  the  quality  of  the  particular  indicators  near 
the  singular  point)  we  also  computed  the  error-indicators  in  the  key  elements  using 
the  extraction  method  as  outlined  above.  In  Table  4a  (resp.  Table  4b)  we  give  the 

results  for  linear  and  quadratic  triangles  for  a  =  -  for  the  mesh  shown  in  Fig.  5a 
^  4 

(resp.  a  =  -  for  the  mesh  shown  in  Fig.  4a).  In  Tables  5a  (resp.  Table  5b)  we 

give  the  results  for  bilinear  (resp.  biquadratic)  squares  for  o  =  ^  and  in  Table  5c 

1 

we  give  the  results  for  a  =  -  and  biquadratic  squaxes.  In  the  above  computations 

4 

for  the  meshes  of  triangles  we  let 


:=  (ij  —  3h,  -I-  3h)  x  (ij  —  3/i,  +  3h)  (52) 

while  for  the  meshes  of  squares  we  used 

'^h  ■—  ~  2^’  ^  (®2  ~  +  3/i)  (53) 

We  observe  that  the  and  estimates  of  the  pollution  have  effectivity 

indices  between  0.9- 1.4. 

We  also  computed  the  A^^^^-estimates  of  the  pollution  for  the  cases  given  in 
Tables  5a,  5b  and  5c.  We  obtained  the  following  effectivity  indices 


=  4.76, 

=  4.83, 

=  4.81 

(54) 

for  the  case  in  Table  5a, 

=  5.68, 

=  4.80, 

=  5.32 

(55) 

for  the  case  in  Table  5b,  and 

=  16.53, 

nf  =  13.43, 

^(3)  ^  14  77 

(56) 
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for  the  case  in  Table  5c.  From  these  results  we  see  that  the  Af^^^-estimate  grossly 
overestimates  the  pollution-error  and  should  not  be  used. 

In  the  Remcirk  3.2  above  we  noted  that  the  derivations  of  the  estimates  for  the 
pollution  are  valid  only  when  the  ERp  error- indicators  are  used  but  nevertheless 
they  could  be  computed  from  other  error  indicators  too  provided  that  these  indi¬ 
cators  are  robust  modulo  the  pollution.  In  Table  3  we  gave  examples  of  pollution 
estimates  based  on  the  ERpB  error-indicators.  To  illustrate  further  the  point  of 
Remark  3.2  we  computed  the  Af  ^^^-estimates  for  the  case  in  Tables  5a,  5b  and  5c 
using  the  ZZ  error-indicators.  We  obtained  the  following  effectivity  indices 


=  1.16, 

= 

1.33, 

=  1.24 

(57) 

for  the  case  given  in  Table  5a, 

/cP  =  2.27, 

2.35, 

=  2.31 

(58) 

for  the  case  given  in  Table  5b,  and 

=  1.78, 

= 

1.68, 

=  1.72 

(59) 

for  the  case  given  in  Table  5c.  As  in  the  case  of  the  ERpB-based  pollution  estimate 
we  get  (reasonable)  overestimation  of  the  error  which  is  related  to  overestimation 
by  the  ZZ  error-indicators  in  the  key  elements. 

5  Adaptive  control  of  the  pollution-error 

Let  t%  be  a  given  tolerance  and  let  us  assume  that  the  goal  is  to  construct  a 
mesh  by  employing  regular  subdivision  of  the  elements  of  an  initial  mesh  T°, 
such  that 


ll|v?'*IIL.<  «ll|vf‘IIL.  (60) 

where  C  T°.  The  mesh  Tf^,  which  achieves  local  control  of  the  pollution-error 
in  according  to  (60)  may  be  obtained  using  the  following  algorithm: 

1.  Let  Tf^  =  T°  and  go  to  3, 

2.  For  each  element  r  £  T/^  do: 

2.1.  Compute  t  £  T^,  t  ^ 

2.2.  If  fi^  >  7max/i^,  subdivide  t. 

3.  Compute  the  finite  element  solution  on  T/^  and  . 
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4.  Check  if 


(61) 


A(,.  < 

If  not,  go  to  2,  otherwise  stop. 

Here  /i,.  is  the  pollution-indicator  for  the  element  t,  Ad is  an  estimator  for 

5  7  ^  (04)  is  a  user-specified  constant.  The  algorithm  may  be  also 
modified  so  that  step  3  is  executed  only  after  a  certain  number  of  refinement  steps. 
In  the  examples  below  we  let  UJ^^  coincide  with  an  element  f  and  we  employed  the 
estimate 


and  for  the  refinement  we  used  the  pollution-indicators 


f^r  =  \J 


(62) 


(63) 


where  x  is  taken  to  coincide  with  a  vertex  of  the  element  f. 

We  employed  the  above  algorithm  with  =  10%  for  meshes  of  bilinear  and 

biquadratic  squares  in  the  X-shaped  domain  (q  =  for  u:^^  ~  f  =  (0.5, 0.625)^ 
(the  element  shown  shaded  gray  in  Figs.  7a,  7b).  We  employed,  as  the  initial 
mesh,  a  uniform  mesh  of  squares  with  Hq  =  -  in  the  second  and  third  quadrants 

1 

and  elements  of  size  fig  =  -  in  =  [0.25,0.75]*.  In  Fig.  7a  (resp.  Fig.  7b)  we 

_  8 

show  the  final  mesh  of  bilinear  (resp.  biquadratic)  squares  constructed  by  the 
algorithm  for  which  we  have  — -  =  9.57%  <  10%  { resp.  — -  =  9.81%  <  10% ). 

In  both  cases  we  employed  the  A4^^^-estimate  of  the  pollution  computed  using  the 
error-indicators  obtained  from  the  extraction  method  in  the  key  elements  and  the 
ERpPp-|-3  error- indicators  in  the  rest  of  the  mesh.  We  also  note  that  in  the  initial 
uniform  mesh  of  bilinear  (resp.  biquadratic)  squares  we  have  \ 

=  0.994  fresp.  t 


resp. 


Vf 


=  0.017 )  while  in  the  final  mesh 


Vf 


=  0.035 
Vf 


Ml  ftIMT  1M  f»IMT  III  fill!  7 

0.923^ .  Thus  by  controlling  the  pollution-error  to  less  than  10%  the  value  of  the 

indicator  we  get  a  very  reliable  local  indicator. 

We  also  considered  the  domaiin  ABCDEFGH  shown  in  Fig.  8  and  exact  solution 


y)  =  rc^  sin  -|-  sin 


(64) 
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where  tc  =  J{x^  -  (xc)i)^  +  (2:2  “  {^0)2)^  =  tan  ( - and  simi- 

VXj  -  (.XpJi/ 

2  1 

larly  for  and  i?£;,  «c  =  o  and  =  -•  We  let  r  be  the  element  shown  shaded  in 

00 

Fig.  8  and  we  adapted  the  mesh  to  achieve  HIVj^Hlf  <  20%|||Vi‘^'‘|||^.  Byu}^^  we  mean 
the  patch  which  includes  f  and  two  mesh- layers  around  it.  We  employed  quadratic 
element  (p  =  2)  and  we  based  the  estimate  M.f  and  the  pollution-indicators 
on  the  ERpB.  In  Figs.  8a  and  8b  we  give  the  pollution-adaptive  meshes  (obtained 
used  7  =  0.8)  for  two  locations  of  the  element  f,  (shown  shaded  gray  with  a  bullet 
at  the  centroid).  In  both  the  examples,  is  the  patch  which  includes  all  the 
shaded  elements  and  one  layer  of  elements  around  them.  In  Fig.  8a  we  show  the 

A4-  77- 

final  grid  for  which  we  have  — -  =  18.6%  <  20%.  For  this  mesh  .  =  0.92 


while  in  the  initial  mesh 


Vf 


a  0.  In  Fig.  8b,  for  the  second  location  of  the 

IlieAlllf 

_ 

element  t,  we  show  the  final  grid  for  which  — -  =  29.8%  <  30%.  For  this  mesh, 

Vf 

0.  We  also  computed  the  effectivity-index  for  the  pollution-estimate 


Vr 


for  the  final  mesh  shown  in  Fig.  8a;  its  value  was  found  1.89.  Hence  we  note 
that  the  pollution-estimates  were  based  on  ERpB  produce  slight  overrefinement  of 
the  mesh. 


6  Summary  of  conclusions 

In  this  work  we  constructed  a-posteriori  estimates  for  the  pollution-error.  The 
main  tools  are  the  local  error-indicators  and  a  finite-element  approximation  of  the 
Green’s  function.  We  also  proposed  an  algorithm  for  the  adaptive  control  of  the 
pollution-error  in  any  patch  of  interest. 

The  major  conclusions  of  this  study  are: 

1.  The  quality  of  the  estimates  of  the  pollution-error  depends  on  the  quality 
of  the  error-indicators  near  the  singular  points.  Existing  estimators  over-  or 
underestimate  the  error  in  the  key  elements  depending  on  the  unsmoothness 
of  the  solution  and  the  degree  of  the  elements.  Hence  it  is  necessary  to 
construct  estimators  with  element  effectivity  indices  close  to  one  near  the 
singular  points.  This  will  be  addressed  in  another  work. 

2.  It  is  possible  to  control  the  pollution-error  in  any  patch  of  elements  by  refining 
the  elements  with  the  biggest  pollution-indicators. 

3.  When  a  direct  solver  is  used  the  cost  of  the  pollution-estimates  is  negligible. 
It  is  practically  the  same  cost  as  computing  with  a  few  additional  right-hand 
sides. 
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4.  For  elements  of  degree  p  >  2  it  is  also  possible  to  construct  estimates  of  the 
pollution  in  the  value  of  the  finite  element  solution.  This  will  be  addressed 
in  a  future  paper. 
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Appendix 

We  will  now  outline  the  recipe  for  the  local  construction  of  equilibrated  element- 
residuals  which  is  employed  in  the  numerical  examples.  This  construction  was 
proposed  in  [21]- [26].  We  let  be  a  function  defined  piecewise  as  a  ^-degree 
polynomial  on  each  edge  e.  We  define  (see  Fig.  2) 


=  r  ,  0‘,  =  -r  (A.l) 

out  in 

Given  any  interior  element  t  and  g,  0  <  q  <  p,  the  aim  is  to  determine  edgewise 
g-degree  polynomials  0^,  such  that 


+■£  I  Kv  =  o  V  t.  €  sJ(t) 


(A.2) 


The  residuals  which  satisfy  (A.2)  are  said  to  be  q-order  equilibrated.  Eq.  (A.2)  is 
equivalent  to 

^  f  0‘r<l>,  == l  =  (A.3) 

sCdr 


where  is  the  t-th  shape-function  of  the  element  and  denotes  the  dimension 
of  5J(r). 

A.l  Equilibration  for  an  interior  vertex  X 

A  set  of  edgewise  polynomial  corrections  for  edges  connected  to  interior  vertices 
is  constructed  as  follows; 

a.  Linear  corrections. 

Let  us  first  determine  edgewise  linear  corrections  in  the  form 


r  =  -1-  r’Va) 


(A.4) 


r,  :=  j^(2A;  -  :=  jl(2A;  -  X[)  (A.5) 

where  A:  =  1,  2  are  the  linear  shape-functions  defined  over  the  edge  e.  Note 
that 

0‘’>’ =  J^0‘XI  ,  it  =1,2  (A.6) 

Let  X  denote  an  interior  vertex  of  the  mesh  with  the  element  and  the  edges 

,  i  =  1,...,4,  connected  to  it,  ^ls  shown  in  Fig.  9a  and  let  denote  the 

local  number  (1  or  2)  of  the  vertex  X  with  respect  to  the  edge  .  The  values  of 
0**  ’"f®* )  are  obtained  from  the  linear  system 
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E  I  Kfi'x  =  . 

-ft  *  * 


fc= 


tCdr^ 


which  reads 


•  1 

0 

0 

-r 

-1 

1 

0 

0 

< 

.  =  < 

0 

-1 

1 

0 

0^3  Me 3) 

.  0 

0 

-1 

1. 

■  -^rxi<f>x)  - 

(A.7) 


(A.8) 


Here  we  cissumed  that  the  directions  of  edge- normals  are  the  ones  shown  in  Fig.  9a. 

The  matrix  has  exactly  one  zero  eigenvalue  with  corresponding  eigenvector 
[1, 1, 1, 1]^.  Moreover  from  the  orthogonality  condition  we  have 


=  (A-9) 

k=l 

and  the  system  is  consistent.  Particular  solutions  can  be  obtained  from  various 
choices  of  the  free-constant.  For  example  one  may  choose  )  =  0  (e.g.  [24]) 

or  one  may  choose  the  constant  by  minimizing  various  norms  of  the  solution  of 
(A.8),  namely 


(A.10) 

k=l 

where  W/,  denotes  the  weight  associated  with  e^.  Here  we  employed  the  choice 
“'fc  =  1^  [25])- 

b.  Higher-order  corrections. 

Let  us  denote  by  k  =  1,...,(9-|-1)  the  basis-functions  which  do  not  vanish 
on  the  edge  e.  We  extend  the  conjugate  basis  introduced  in  (A. 5)  as  follows: 

r  0,  if  j  <  i 

/  =  i  (A.ll) 

ll,  ifi  =  l 

where  3  <  i,j  <  (?  +  1).  After  the  linear  edge  corrections  have  been  determined, 
the  higher-order  corrections  can  be  computed  from 


z  =  3,...,(g-t-l),  sCdr^t  (A.12) 
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Note  that  the  higher-order  corrections  can  be  determined  explicitly  and  are  defined 
uniquely  for  each  edge. 

A. 2  Equilibration  for  a  vertex  X  on  the  boundary  of  the  domain. 

The  edgewise  polynomial  corrections  for  edges  connected  to  boundary-nodes 
are  constructed  as  follows; 

a.  Linear  corrections. 

Let  X  be  a  vertex  on  the  boundary  of  the  domain,  with  elements  f  =  1, . . . ,  3 
and  edges  ef ,  j  =  1, . . . ,  4  connected  to  it,  as  shown  in  Fig.  9b.  The  edges  and 
lie  on  the  boundary.  For  the  vertex  X  we  get  the  following  linear  system 


Y. 


which  reads  as. 


■-1  1  0 

0  -1  1 

.  0  0-1 


•  0^1  ' 

Qe^M‘2) 

•< 

J 

k  =  1,2,3. 


■  -^rX{4>x)  ‘ 


(A.13) 


(A.14) 


The  directions  of  the  edge  normals  are  shown  in  Fig.  9b.  We  will  have  various 
possibilities  for  the  type  of  boundary  conditions  on  edges  ef  and  Based  on 
these  possibilities  a  solution  to  the  above  system  is  obtained.  The  following  types 
of  boundary  conditions  are  possible: 

(i)  Dirichlet  boundary  conditions  on  both  edges  and 

The  system  (A.14)  has  an  infinite  set  of  solutions.  A  particular  solution 
which  minimizes  (A.  10)  is  selected. 

(ii)  Neumann  boundary  conditions  on  both  edges  ef  and  e^: 

The  system  (A.14)  has  a  unique  solution  which  is  obtained  by  setting  ^ 

aud  equal  to  zero. 

(iii)  One  of  the  edges  is  on  the  Dirichlet  boundary  while  the  other  is  on 

the  Neumann  boundary: 

A  unique  solution  is  obtained  for  the  system  (A.14)  by  setting  equal 

to  zero  for  the  edge  on  the  Neumann  boundary. 

b.  Higher-order  corrections. 

The  higher-order  corrections  are  computed  following  exactly  the  same  steps  as 
for  the  interior  vertex. 
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List  of  Figures 

Fig.  1.  An  example  of  polygonal  domain. 

Fig.  2.  An  interior  edge  with  its  unit-normal  and  the  elements  and  attached 
to  it. 

Fig.  3.  An  example  of  a  patch  of  elements  included  in  a  larger  patch  which 
includes  and  three  mesh-layers  around  it. 

Fig.  4.  (a)  The  i-shaped  domain  covered  by  a  mesh  of  triangles  with  ^  ~  Jg! 

(b)  Detail  of  the  mesh  showing  the  enumeration  of  the  elements  in  the  neighborhood 
of  the  singular  point. 

Fig.  5.  (a)  The  slit-domain  covered  by  a  three-directional  mesh  of  triangles  with 

h  =  — ;  (b)  Detciil  of  the  mesh  near  the  singular  point  showing  the  enumeration 
16 

of  the  elements. 

Fig.  6.  (a)  The  X-shaped  domain  covered  by  a  mesh  of  squares  with  ^ 

(b)  Detail  of  the  mesh  in  the  neighborhood  of  the  singular  point  showing  the 
enumeration  of  the  elements. 

Fig.  7.  Adaptive  control  of  the  pollution  error;  Final  mesh  obtained  for  the  control 
of  the  pollution-error  in  the  shaded  element  f  for  t%  =  10%.  X-shaped  domain 

with  a  =  -.  Adaptive  mesh  of: 

M- 

(a)  Bilinear  squares  (p  =  1)  which  satisfies  — -  =  9.57%  <  10%. 

Vr 

M- 

(b)  Biquadratic  squares  (p  =  2)  which  satisfies  — -  =  9.81%  <  10%. 

Vf 

Fig.  8.  Adaptive  control  of  the  pollution  error:  Adaptive  mesh  of  quadratic  tri¬ 
angles  (p  =  2)  obtained  for  the  control  of  the  pollution-error  in  the  shaded  ele- 

ment,  f  with  a  bullet  at  the  centroid.  Polygonal  dommn,  u{x,y)  =  (3^0)  + 

4  sin 


(a)  Adaptive  mesh  which  satisfies  — -  =  18.6%  <  20%; 

Vf 

M 

(b)  Adaptive  mesh  which  satisfies  — ~  =  29.8%  <  30%. 

Vf 

Fig.  9.  (a)  An  interior  vertex  X,  (b)  A  boundary  vertex  X.  Here  we  show  the 
elements  and  the  edges  and  the  unit-normals  assigned  to  the 

edges. 
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Effectivity  indices  in  the  elements  in 
the  neighborhood  of  the  singular  point 

Uniform  mesh  of  linear  triangles  (p  = 

1) 

Element 

lillll 

ERpPp  + 1 

ERpB 

ZZ 

PRpPp-¥3 

No. 

mg 

u(r,d)  =  risin(^) 

iH 

2.544 

1.897 

0.475 

0.759 

3.008 

1.959 

0.562 

0.788 

Bfil 

2.737 

1.656 

0.502 

0.916 

0.106 

1.930 

1.361 

0.603 

1.031 

5 

0.138 

1.787 

1.037 

0.701 

1.098 

6 

0.114 

0.729 

0.558 

1.298 

1.178 

u(r,d)  =  risin(Y) 

1 

0.129 

1.696 

1.268 

1.028 

0.964 

2 

1.372 

1.037 

1.183 

0.979 

3 

0.127 

1.306 

0.776 

0.907 

1.026 

4 

0.127 

1.585 

0.920 

0.907 

1.026 

5 

1.739 

1.269 

1.183 

0.979 

6 

0.129 

1.972 

1.473 

1.028 

0.964 

Table  la.  The  quality  of  the  error  indicators  in  the  neighborhood  of  the  singular 
point:  The  element  effectivity  indices  for  the  mesh  and  elements  shown  in  Figs.  4a 

and  4b,  respectively.  Uniform  mesh  of  linear  trieuigles  (p  =  1),  A  =  — ,  u(r,r9)  = 


r“sin(ai?)  (a  =  ^)- 


16’ 
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Effectivity  indices  in  the  elements  in 
the  neighborhood  of  the  singular  point 

Uniform  mesh  of  linear  triangles  (p  = 

1) 

Element 

No. 

miB 

ERpPp+  1 

ERpB 

ZZ 

PRpPp  +  3 

iny 

ii(r,t9)  =  risin(^) 

1 

0.111 

2.861 

2.133 

2 

0.099 

3.498 

2.250 

3 

0.190 

2.071 

1.667 

WmM 

4 

0.108 

2.625 

1.575 

5 

0.079 

1.718 

1.231 

0.868 

6 

0.172 

1.076 

0.687 

1.009 

■1 

«(r,j?)  =  r 

isin(|) 

1 

0.106 

2.003 

1.494 

0.775 

0.886 

2 

2.074 

1.432 

1.169 

0.929 

3 

0.187 

1.266 

0.799 

1.154 

1.035 

4 

0.907 

0.724 

1.133 

1.077 

5 

0.113 

2.060 

1.200 

1.120 

1.010 

6 

0.174 

1.221 

1.181 

0.933 

0.898 

Table  lb.  The  quality  of  the  error  indicators  in  the  neighborhood  of  the  singular 
point;  The  element  effectivity  indices  for  the  mesh  and  elements  shown  in  Figs.  5a 

and  5b,  respectively.  Uniform  mesh  of  linear  triangles  (p  =  1),  ,  u(r,  i?)  = 

11 

r"sin(at?) 
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Effectivity  indices  in  the  elements  in 
the  neighborhood  of  the  singular  point 

Uniform  mesh  of  quadratic  triangles  (p 

=  2) 

Element 

No. 

mill 

ERpPp  + 1 

ERpB 

ZZ 

PRpPp  +  3 

u(r,d)  =  rism(^) 

2.370 

1.915 

0.513 

0.593 

2.655 

2.153 

0.574 

0.657 

2.397 

1.914 

0.572 

0.779 

0.147 

1.776 

1.441 

0.662 

0.898 

5 

0.155 

1.407 

1.086 

0.848 

0.986 

6 

0.146 

0.681 

0.493 

1.266 

1.046 

u(r,i?)=  rJsin(Y) 

1 

1.543 

1.275 

1.178 

0.833 

2 

1.361 

1.084 

1.256 

0.939 

3 

1.002 

0.713 

1.079 

1.037 

4 

0.171 

1.201 

0.881 

1.079 

1.037 

5 

0.148 

1.650 

1.325 

1.256 

0.939 

6 

0.170 

1.791 

1.474 

1.178 

0.833 

Table  Ic.  The  quality  of  the  error  indicators  in  the  neighborhood  of  the  singular 
point:  The  element  effectivity  indices  for  the  mesh  and  the  elements  shown  in 

Figs.  4a  and  4b,  respectively.  Uniform  mesh  of  quadratic  elements  (p  =  2),  /i  = 

1  12 

— ,  u(r,t9)  =  r®sin(ai?)  (a  =3^3)- 
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Effectivity  indices  in  the  elements  in 
the  neighborhood  of  the  singular  point 

Uniform  mesh  of  quadratic  triangles  {p 

=  2) 

Element 

No. 

ERpPp  + 1 

ERpB 

ZZ 

PRpPp  +  3 

IIKIII^ 

u(r,  d)  =  r*  sin(^) 

1 

0.123 

2.689 

2.166 

0.336 

0.456 

2 

0.119 

3.116 

2.521 

0.471 

0.509 

3 

0.222 

1.851 

1.743 

0.588 

0.665 

4 

0.133 

2.291 

1.829 

0.703 

0.823 

5 

0.103 

1.610 

1.301 

0.943 

0.924 

6 

0.197 

0.945 

0.679 

1.108 

1.022 

«(r,t?)  =  r 

i  sin(^) 

1 

0.125 

1.820 

1.485 

1.895 

0.732 

2 

0.114 

1.824 

1.471 

msm 

0.828 

3 

0.234 

1.130 

0.798 

1.009 

4 

0.117 

0.866 

0.640 

1.046 

5 

0.128 

1.698 

1.326 

1.289 

0.957 

6 

0.232 

1.044 

1.019 

0.775 

Table  Id.  The  quality  of  the  error  indicators  in  the  neighborhood  of  the  singul£ir 
point:  The  element  effectivity  indices  for  the  mesh  and  the  elements  shown  in 
Figs.  5a  and  5b,  respectively.  Uniform  mesh  of  quadratic  elements  {jp  =  2),  h  = 

— ,  u(r,i?)  =  r“sin(aT?)  (q  =  2)’ 
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EfFectivity  indices  in  the  elements  in  the  neighborhood  of  the  singular  point 


Uniform  mesh  of  bilinear  squares  (p  =  1) 

u(r,  t?)  = 

risin(|) 

Element 

No. 

■iiBniiai 

mmm 

ERpPp  +  1 

ERpPp  +  2 

ERpPp  +  3 

ERpPp  +  4 

PRpPp  +  3 

ZZ 

1 

0.005 

0.667 

0.709 

0.796 

0.807 

0.811 

0.408 

2 

0.008 

0.325 

0.471 

0.487 

0.509 

0.510 

0.403 

3 

0.008 

0.821 

0.864 

0.923 

0.933 

0.892 

0.686 

4 

0.005 

0.527 

0.538 

0.557 

0.559 

0.600 

0.477 

5 

0.006 

1.240 

1.358 

1.443 

1.468 

1.219 

1.107 

6 

0.287 

0.770 

0.859 

0.990 

1.012 

1.106 

0.795 

7 

0.313 

0.798 

0.845 

1.008 

1.020 

1.014 

0.727 

8 

0.009 

2.389 

2.456 

2.904 

2.922 

2.716 

0.627 

9 

0.007 

0.883 

1.057 

1.110 

1.141 

1.207 

0.887 

10 

0.262 

1.084 

1.212 

1.402 

1.434 

1.381 

1.555 

11 

0.006 

0.757 

0.816 

0.938 

0.954 

0.974 

0.347 

12 

0.006 

1.284 

1.457 

1.540 

1.574 

1.613 

2.379 

Table  2.  The  quality  of  the  error  indicators  in  the  neighborhood  of  the  singular 
point:  The  element  eflFectivity  indices  for  the  mesh  and  the  elements  shown  in 

Fig.  6a  and  6b,  respectively.  Uniform  mesh  of  bilinear  squares  (p  =  1),  ft  =  — , 
1 

«(r,i?)  =  rs  sin(|). 
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A-posteriori  estimates  of  the  pollution-error 
Estimates  based  on  the  ERpB  estimator 
Mesh-patch  centered  at  z  =  (0.5, 0.5) 
u(r,»>)  =  risin(|) 

Exact 

Quantity  v - - - 

V  V  /c(€)  v  «(£!) 


Uniform  mesh  of  linear  triangles  (p  =  1) 


.0332 

.0631 

1.90 

.0640 

1.92 

.0281 

.0585 

2.08 

.0597 

2.13 

.0435 

.0860 

1.98 

.0875 

2.01 

Uniform  mesh  of  linear  triangles  (p  =  2) 

.0157 

.0382 

2.43 

.0382 

2.43 

.0132 

.0349 

2.64 

.0349 

2.64 

.0205 

.0517 

2.52 

.0517 

2.52 

Table  3 

exact  values  of 


The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The 


dxi 


(*)|,  I  ■(«)[  or  IWj'*"'*!!*),  the  corresponding  and 


Af^^^-estimates  (computed  by  employing  the  error-indicators  from  the  ERpB  esti¬ 
mator)  and  the  corresponding  effectivity  indices.  Uniform  mesh  of  linear  (p  =  1) 

or  quadratic  {p  =  2)  triaingles  shown  in  Fig.  4a  [h  =  — u(r,i?)  =  sin(— ). 

'  In/  '  51 
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A-posteriori  estimates  of  the  pollution-error 


Estimates  based  on  the  extraction  method 


Mesh-patch  Cj|^  centered  at  z  =  (0.5, 0.5) 


u(r,d)  =  risin(^) 


Exact 


V  V  «(€)  V  k(v) 


Uniform  mesh  of  linear  triangles  (p  =  1) 


Quantity  v 


.0498 

.0459 

0.92 

.0527 

.0437 

.0441 

1.01 

.0481 

.0663 

.0637 

0.96 

.0714 

& 

|W^|(i) 


Uniform  mesh  of  quadratic  triangles  (p  =  2) 


It'*’! 

.0293 

.0276 

0.94 

.0307 

1.05 

.0258 

.0252 

0.98 

.0265 

1.03 

.0390 

.0374 

0.96 

.0406 

1.04 

Table  4a.  The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The 

exact  values  of  — (a)  ,  -r-— (ic)  or  |W^'*|(x),  the  corresponding  and 

'  C/X-y  *  '  0X2 

A^^^hestimates  (computed  by  employing  the  error-indicators  from  the  extraction 
method  in  the  key  elements)  and  the  corresponding  effectivity  indices.  Uniform 

mesh  of  triangles  shown  in  Fig.  5a,  h  =  — ,  u(r,  1?)  =  ri  sin(— ).  The  error- 

lo  4 

indicators  computed  from  the  extraction  method  were  employed  in  one  mesh-layer 
around  the  singular  point;  the  error-indicators  from  ERpPp  -|-  3  were  employed  in 
the  rest  of  the  mesh. 
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A-posteriori  estimates  of  the  pollution-error 


Estimates  based  on  the  extraction  method 


Mesh-patch  centered  at  z  =  (0.5, 0.5) 


u(r,d)  =  ri8m(|) 


Quantity  v 


Exact 


k(v) 


K(t5) 


Uniform  mesh  of  linear  triangles  (p  =  1) 


.0332 

.0343 

— 

1.03 

.0382 

1.15 

.0281 

.0338 

1.18 

.0379 

1.35 

|W/‘1(Z) 

.0435 

.0482 

1.11 

.0538 

1.24 

Uniform  mesh  of  quadratic  triangles  (p  =  2) 

^<*'1 

.0157 

.0152 

0.97 

.0161 

1.03 

C'*'l 

.0132 

.0149 

1.13 

.0149 

1.13 

IW“‘|(i) 

.0205 

.0209 

1.02 

.0219 

1.07 

Table  4b.  The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The  ex¬ 


act  values  of 


lavC' 


(aj)j,  I  ^  -(g)!  or  |W^'‘|(a5)  and  the  corresponding  and 

'  *  vX2 

A4(*)-e8timates  (computed  by  employing  the  error-indicators  from  the  extraction 
method  in  the  key  elements)  and  the  corresponding  effectivity  indices.  Uniform 

mesh  of  triangles  shown  in  Fig.  4a,  h  =  — ,  «(r,  t?)  =  r?  sin(— ).  The  error- 

10  ^  O 

indicators  computed  from  the  extraction  method  were  employed  in  one  mesh-layer 
arotmd  the  singular  point;  the  error-indicators  from  ERpPp  -I-  3  were  employed  in 
the  rest  of  the  mesh. 
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A-posteriori  estimates  of  the  pollution-error 


Estimates  based  on  the  extraction  method 


Uniform  mesh  of  bilinear  squares  (p  =  1) 


u(r,t?)  =  risin(|) 


Mesh-patch  centered  at  x  =  (0.5, 0.5) 


Quantity  Vj 

Exact 

V 

k(i)i) 

k(vi) 

V-(*) 

0.0212 

0.0242 

1.14 

0.0301 

1.41 

dll  1 

dVf'*  1 

1  * 

0.0188 

0.0199 

1.06 

0.0228 

1.21 

dxj  1 

0.0288 

0.0314 

1.09 

0.0378 

1.31 

Table  5a.  The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The 


exact  values  of  [  —  (g)[,  ■■(»)[  or  [VV^^Kx),  the  corresponding  and 


dx. 


dx 


A^^^^-estimates  (computed  by  employing  the  error- indicators  from  the  extraction 
method  in  the  key  elements)  and  the  corresponding  elfectivity  indices.  Uniform 

mesh  of  bilinear  squares  (p  =  1),  h  =  — ;  u(r,i?)  =  rs  sin(— ).  The  error-indicators 

lo  3 

computed  from  the  extraction  method  were  employed  in  two  mesh-layers  around 
the  singular  point;  the  error-indicators  from  ERpPp  -f  3  were  employed  in  the  rest 
of  the  mesh. 
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A-posteriori  estimates  of  the  pollution-error 


Estimates  based  on  the  extraction  method 


Uniform  mesh  of  biquadratic  squares  (p  =  2) 


u(r,  d)  =  risin(|) 


Mesh-patch  centered  at  z  =  (0.5, 0.5) 


Quantity  Uj 

Exact 

V 

»'i 

«(Vi) 

K(t'i) 

0.0168 

0.0202 

1.20 

0.0233 

1.38 

0.0154 

0.0177 

1.15 

0.0184 

1.19 

|Vvf‘|(*) 

0.0228 

0.0269 

1.18 

0.0297 

1.30 

Table  5b.  The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The 
exact  values  of  [  ^  (g)[,  [  (x)[  or  the  corresponding  and 


A^^^)-estimates  (computed  by  employing  the  error-indicators  from  the  extraction 
method  in  the  key  elements)  and  the  corresponding  effectivity  indices.  Uniform 


1  i  .  ” 

mesh  of  biquadratic  squares  (p  =  2),  A  =  u(r,  i?)  =  rs  sin(— ).  The  error- 

o  3 

indicators  computed  from  the  extraction  method  were  employed  in  two  mesh-layers 
Mound  the  singular  point;  the  error-indicators  from  ERpPp  4-  3  were  employed  in 
the  rest  of  the  mesh. 
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A'posteriori  estimates  of  the  pollution-error 


Estimates  based  on  the  extraction  method 


Uniform  mesh  of  biquadratic  squares  (p  =  2) 


u(r,d)=  rTsin(|) 


Mesh-patch  centered  at  z  =  (0.5, 0.5) 


^ _ 

Exact 

V 

*'1 

0.0272 

0.0332 

1.22 

0.0429 

1.58 

0.0323 

0.0397 

1.23 

0.0485 

1.51 

0.0420 

0.0517 

1.23 

0.0647 

1.54 

Table  5c.  The  quality  of  the  a-posteriori  estimates  of  the  pollution-error:  The 
exact  veJues  of  (g)|  or  the  corresponding  and 


A^^^^-estimates  (computed  by  employing  the  error-indicators  from  the  extraction 
method  in  the  key  elements)  and  the  corresponding  elfectivity  indices.  Uniform 


mesh  of  bilinear  squares  (p  =  2),  A  =  u(r,  i?)  =  r*  sin(— ).  The  error-indicators 

computed  from  the  extraction  method  were  employed  in  two  mesh-layers  around 
the  singular  point;  the  error-indicators  from  ERpPp  -f  3  were  employed  in  the  rest 
of  the  mesh. 
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Fig.  4a 


Fig.  4b 
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The  Laboratory  for  Numerical  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Tedmology  of  the  University  of  Mar^and,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals: 

To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  empha^  on  the  numeric^  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  govenunent  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

To  assist  with  the  education  of  numerical  analysts,  espedally  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisdplinaiy  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Sdence  Departments.  This  includes  active 
collaboration  with  government  agencies  sudi  as  the  National  Institute  of  Standards  and 
Tedmology. 

To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agendes 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  L  Babujka,  Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Sdence  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


